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The current standard model of cosmology, ACDM, requires dark matter to make up around 25% 
of the total energy budget of the Universe. Yet, quite puzzlingly, there appears to be no candidate 
particle in the current Standard Model of particle physics. Assuming the validity of the cold dark 
matter (CDM) paradigm, dark matter has evaded detection thus far either because it is intrinsically 
a weakly interacting substance or because its interactions are suppressed by its high constituent 
mass and low number density. Most approaches to explain dark matter to date assume the former 
and therefore require beyond-the-Standard-Model particles that have yet to be observed directly or 
indirectly. Given the dearth of evidence for this class of candidates it is timely to consider the latter 
possibility, which allows for candidates that may or may not arise from the Standard Model. In this 
work we extend a recent study of this general class of so-called macro dark matter-candidates with 
characteristic masses of grams and geometric cross sections of cm^. We consider new bounds that 
can be set using existing data from the resonant bar gravitational wave detectors NAUTILUS and 
EXPLORER. 


I. INTRODUCTION 

The last two decades have ushered in a new era of pre¬ 
cision cosmology with a plethora of modern experiments 
and observations all leading to the so-called ACDM con¬ 
cordance model of the Universe. A crucial component 
of this model is the cold dark matter (CDM) making up 
around 25% of the total energy density of our Universe, 
as inferred from observations on large scales (see e.g. Hi¬ 
ll])- On intermediate scales its existence is observed indi¬ 
rectly due to the discrepancies between the gravitational 
and luminous masses of large astrophysical systems. The 
dark matter cannot be made of ordinary baryonic matter 
(see e.g. [SHZj, and references therein) nor can it consist 
of any fundamental particle of the Standard Model of 
particle physics. 

Although a tremendous effort has been devoted to the 
study of particle candidates that are intrinsically weakly 
interacting, such as weakly interacting massive particles 
and axions, alternative dark matter candidates deserve 
increased attention given the lack of direct detection of 
those canonical possibilities to date (see e.g. 0 E])- 
Other candidates can couple to Standard Model parti¬ 
cles with high probability but are nevertheless effectively 
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weakly interacting or “dark” because they are very mas¬ 
sive, and therefore have a lower number density. Several 
examples exist in the literature, most notably nuclear- 
dense candidates with a Standard Model basis (e.g. uni- 
ill]). Other similarly massive candidates include compact 
objects with some connection to the Standard Model (e.g. 
[IMS]), primordial black holes |16j . and other candi¬ 
dates that may be found in the literature (e.g. dZHU]). 

A wide range of Earth-based constraints on nuclear- 
dense candidates (nuclearites) was considered in [22] and 
more recent bounds were presented in [H]. Such con¬ 
straints are only made for candidates with a fixed inter¬ 
nal mass density of roughly 3.5 x 10^"^ g/cm^, and there¬ 
fore they obey the specific cross section-mass relation 
CTx cx Recently, however, a model-independent 

study of the general class of “macroscopic” candidates 
has been considered in |24] . These macro dark mat¬ 
ter candidates, referred to as Macros, have characteristic 
masses and geometric cross sections of grams and cm^ 
and may or may not have a Standard Model origin. A 
number of Earth-based, astrophysical, and cosmological 
observations were used in |24j to place constraints on 
Macros. 

Here we improve upon the constraints found in |24j by 
extending existing work on nuclearite constraints from 
resonant bar gravitational wave detectors for applica¬ 
tion to macro dark matter. The use of large aluminum 
bars as detectors of gravitational waves has been stud¬ 
ied for more than five decades [25] . and their sensitivity 
to cosmic rays has been appreciated for some time [26] . 
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Because of the thermoacoustic effect and their resonant 
properties, such aluminum bars are also sensitive (when 
cryogenically cooled) to other exotic cosmic impactors, 
such as monopoles m and nuclearites [28]. The reso¬ 
nant bar experiments NAUTILUS [H], which ran from 
2003-2012 in Frascati, Italy, and EXPLORER [30], which 
ran from 2003-2009 at CERN, were recently used to con¬ 
strain nuclearite dark matter with masses up to ~ 10 “^ 

g in]- 

Compared to nuclearites, macro dark matter has the 
added complexity that both the cross section and mass 
are independent parameters, and both of these determine 
Macro detectability. Macros of a fixed mass with larger 
cross sections, for example, will experience more drag and 
so would arrive at a detector with less energy; they would 
therefore have a lower probability to reach a detection 
threshold. Here we expand on the work of m and |35| . 
to compute the detectability of macro dark matter as a 
function of both its cross section and mass. We then use 
the null results of m to constrain a portion of the Macro 
parameter space and improve on the constraints found in 
[33] on the elastic scattering cross section of macro dark 
matter with baryons. 


II. MACRO DETECTION VIA RESONANT 
BAR DETECTORS 

A. Macro velocity evolution 

Given an incoming Macro galactic velocity, vq, its evo¬ 
lution as it impacts the Earth as a function of the reduced 
cross section, is approximately given by 

i;x(r) = uo (1) 

where (pr) is the column density encountered by the 
Macro. This column density can have three distinct con¬ 
tributions, the first being atmospheric. Simple atmo¬ 
spheric models give an atmospheric density that depends 
exponentially on altitude: 

p{z) = , ( 2 ) 

where z is the altitude above sea level, po — 10 “^ g/cm^, 
and the scale height iL ~ 10 km. Writing z = r cos 9, 
where 9 is the impact angle, the column density encoun¬ 
tered through the atmosphere is 

poo 

(p?’)atm = / dr'p 
J r 

^ p^ 

cos9 ■ ^ ’ 

In addition, we shall conservatively assume that at angles 
9 > it/ 3 the Macro would encounter an additional column 
density of 100 meters water equivalent (MWE), or 

(pr)ob.t = lO^g/cm^ , (4) 


due to surrounding structures, mountains, or other ob¬ 
stacles. Lastly, at angles d > 7 r /2 the Macro would pass 
though the Earth and encounter a column density, 

(pr)® = 2Rep® Icosdl . (5) 


B. Thermo-acoustic detection of cosmic particles 


Upon its passage through the (identical) NAUTILUS 
or EXPLORER detectors, a very dense object would 
deposit energy in a line along its track in the detec¬ 
tor. Because the Macro velocity is supersonic in the alu¬ 
minum, the energy deposition would be nearly instan¬ 
taneous. The resulting near-instantaneous thermal ex¬ 
pansion would source a pressure wave that would excite 
the bar’s longitudinal vibrational modes. As described 
in e.g. [37113^ |3T], for a bar of radius R and length L, 
the energy of the lowest mode is 


where 7 is the Griineisen parameter, p is the bar density. 
Vs is the longitudinal sound speed, and the energy lost 
by a nuclearite or, more generally, a Macro is 


dEy_ 

dx 


payrvl. 


( 7 ) 


The geometric function that depends on the track 
through the bar is 


G{zo, lo, do) 


[ ^ • 

^7rzo^ 

sin cos 

do)] 

ttR 

^ L y 

cos do 



( 8 ) 


where Zq is the distance of the track midpoint to one 
end of the bar, 9o is the angle between the track and 
the bar axis and Iq is the length of the track. For the 
bars described in |3T], i? = 0.3 m and L = 3 m; the full 
geometric acceptance is therefore 19.54 m^ sr. 


C. Application of nuclearite analysis to Macros 


In Ref. |3T] it is reported that, for the parameters 
specific to the NAUTILUS and EXPLORER aluminum 
bars, a vertical impact at the bar center results in an 
energy deposition (in units of temperature) into the fun¬ 
damental mode of 


AE = 10.7 K 


( v.gM) 

\ 10“3c 


4 


( 9 ) 


where, in our notation, ^(M) = (M/1.5ng)^/^ is a func¬ 
tion that characterizes the radius of a passing nuclearite 
as a function of its mass^. 


^ For an object of constant density (Tx oc hence from § 

and 0 the dependence in 0 can be understood. 
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It is straightforward to apply the analysis in [ST] to 
macro dark matter, considering a nuclearite of mass 1.5 
ng would have a cross section (Tq — tt x 10“^® cm^ at 
the reference density of 3.5 x 10^^ g/cm^. The energy 
deposition by a Macro can then be readily translated by 
making the replacement 

C(M)2 ^ ^ (10) 

o’o 

in ([^. We will furthermore use 250 km/s as a reference 
velocity and reinsert the geometric function, G, writing 
the excitation energy as 

Given the detection threshold of 2 K used in [31] and the 
Macro velocity upon impact given in 0 , the fraction of 
impacts for a given pair of Macro parameters (ctx, Mx) 
that exceeds this threshold is determined by performing 
a Monte Carlo simulation over the possible bar impact 
points and Macro trajectories. 

We can, however, give a semianalytic estimate of the 
constrained region in the Ux — Mx plane by recognizing 
that at large a^/My^ the only Macros that can make it 
to the detector with sufficiently high velocity (quanti¬ 
fied below) will come from nearly directly overhead. For 
those impacts G{zq,Iq,9q) is maximized at unity. The 
requirement for resonant bar detection is then 


^^0 


250 km/s 



> 0.62. 


( 12 ) 


This inequality is saturated at the critical value Ux.cj 
which may be solved implicitly as a function of M^: 


CTx.c = - 


Mx 

2{pr) 


W{x). 


(13) 


Here 


X ~ —1.24 X {pr] 


o-Q 


^^0 


Mx \250 km/s 


(14) 


and W {x) is known as the Lambert-W function, defined 
implicitly by the relation 


X = W{x) exp W(x). 


(15) 


W(x) has an infinite number of branches; however, the 
branch W-i{x) delineates the top of the constrained re¬ 
gion in the CTx — -^x plane according to (131. 


We must also ensure that the Macro velocity exceeds 
the sound velocity of aluminum, Vs — 2 x 10“® c; other¬ 
wise the energy loss rate formula ([^ breaks down [31]. 
From 0 this requires that 


c^x < 7 ^ log —Mx 
{pr) Vs 

< 3.7 X 10-3 cm2 f ^ 

VI g 


(16) 


where we have used an approximate atmospheric depth^ 
of (pr) = 103 g/cm^ and vq = 250 km/s. For masses 
larger than roughly 10“ ^^g this turns out to be more 
restrictive than using (131, so it sets the upper edge of 


the region of constraint in that range. 

Lastly, we consider the Macro flux limitation. At 
large CTx/Mx, we cannot use the full geometric accep¬ 
tance of ~ 19.54 m^ sr because of significant additional 
drag that would be experienced by the Macros during 
their passage through (i) the environment surrounding 
the detector (e.g. other buildings, mountains, etc.) and 
(ii) the integrated depth through the Earth. Likewise, 
Macros with small Ux/Mx are more likely to make it to 
the detector with high enough velocity and energy de¬ 
position. Therefore, for the largest a^/My we use only 
1/4 of the geometric acceptance since we can only trust 
their ability to make it through the atmosphere. For 
CTx/Mx < 10“"^ cm^/g a Macro could easily pass through 
surrounding obstacles so that using 1/2 of the full ac¬ 
ceptance is justified. Finally, for Ux/Mx < {p^R<^)~^ ~ 
3 X 10-30 cm^/g the Macro could pass freely through the 
Earth and in that case we may use the full acceptance. 
These approximations ignore the detector efficiency as a 
function of impact trajectory; however this is taken into 
account in the Monte Carlo simulation. 

Given an isotropic dark matter flux of 


—rixUo —1.2 X 10 ®m ^sr ^day ^ 
47r 


1 g 
Mx 


Px 


0.4GeV/cmOy y250 km/s 


■^0 


(17) 


and the combined 3,921 live-time days of NAUTILUS and 
EXPLORER, the maximum number of events expected 
is 


iV=10-(i|) (18) 

for a local dark matter density, px = O.dGeV/cm^, and 
Vq = 250 km/s. A Macro impact would be a random 
(Poisson) process; therefore a null detection of Macro 
passages for which N > 3 events are expected will rule 
them out at greater than 95% confidence. It follows that, 
accounting for the usable fraction of the geometric accep¬ 
tance mentioned above, the following macro dark matter 


^ We have neglected the additional integrated depth the Macro 
must pass through due to concrete in the roofing of the building 
in which EXPLORER was housed, as well as the experimental 
components mounted to the top of the bars m- On the other 
hand, EXPLORER was itself run at 430 m above sea level at 
CERN, so we have overestimated the atmospheric depth to it. 
These two complications (at least) partially offset each other 
and, in any case, they can be neglected considering the level 
of accuracy required here. 
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candidates are ruled out by this analysis: 


[Sx 10-5g, 
Mx < i 2 X 10-5g, 
[l0-5g, 


if ^ < 10-®£S!_. 

Mx ^ g ’ 

if 10-9£S1^ < ^ 

g ^ Mx -- g ’ 

if 10-42JS- < ^ < 10-32J2-. 

g ^ Mx ^ g 


(19) 

Note that at the smallest a^/M^ our bound differs 
slightly from the one found in m because of both our 
choice of vq = 250 km/s and 95% confidence require¬ 
ment. Deviations from these estimates are expected at 
low enough CTx because of the decreasing excitation en¬ 
ergy as a function of CTx, as seen in 0- The results of 
our more accurate Monte Carlo analysis are illustrated 
in Fig. [2 where, however, we have used the fixed value 
of Vo = 250 km/s. 



FIG. 1: Existing constraints (described in detail in [53]) on 
macro dark matter that scatters elastically with baryons and 
our new constraints from resonant bars displayed in red. Note 
the marginal overlap between the resonant bar and large-scale 
structure constraints and that slight improvements to the Sky- 
lab and mica constraints have been made as compared to [53] ; 
both of these points are explained in the Discussion section. 
The CMB constraint [241133] . applicable if Macros couple to 
photons, rules out Macros above the dashed line. 


III. DISCUSSION 


taken in combination with the mica constraints, com¬ 
pletely overlap with these constraints. Nevertheless, our 
results are relevant since they directly constrain Macro- 
baryon coupling which could, in principle, be different 
from Macro-photon coupling. Also, the resonant bar de¬ 
tectors have acted as a local probe that directly con¬ 
strains Macro properties in the solar neighborhood with¬ 
out relying on the (albeit reasonable) assumption that 
the local dark matter is of similar form and composition 
to dark matter in different parts and epochs of the Uni¬ 
verse. 

Here the velocity distribution of Macros was taken to 
be isotropic; in reality the Earth (with the Solar System) 
moves through the Galaxy at approximately 200-250 
km/s and this should result in an anisotropic distribution 
observed on the Earth. Also, the Macro velocity distri¬ 
bution presumably has a tail from which it is common to 
draw an initial Macro velocity, vq significantly larger than 
250 km/s. This means that the constraints inferred here 
are conservative and the region of constraint should ex¬ 
tend upwards to larger Ux, be. the constraint is stronger 
than (Tx/A/x < 3.7 x 10“^ cm^/g. Therefore, greater 
overlap is expected with those determined from large- 
scale structure, which are Ux/Mx < 3.3 x 10“^ cm^/g 
[34] . On the other hand, the range of tJx/Afx constrained 
here is only log-sensitive to Vq, so we do not expect this 
to significantly enhance our results. 

Here we have also made modest improvements to the 
constraints inferred from both Skylab [35] and ancient 
mica [ 35 ] compared to those presented in [ 21 ]. For both 
detectors we have more carefully calculated their accep¬ 
tance, taking into account the column density encoun¬ 
tered by a Macro as a function of its impact angle-this 
resulted in more rounded corners of the constrained re¬ 
gions. For the mica samples studied in [ 35 ] the orienta¬ 
tion during their ~ 500 Myr exposure time is unknown; 
we therefore computed the minimum region that is ruled 
out for a flat mica sample oriented either parallel or 
perpendicular to the Earth’s surface and have presented 
those constraints in Fig. [^ 

The new constraints presented here on macro dark 
matter from resonant bar detectors provide overlap with 
both the Skylab and mica constraints, and also fill in 
part of the gap between the large-scale structure and 
mica bounds on the coupling of Macros to baryons. 
In future studies, however, a dedicated experiment of 
this type is far from ideal due to its relatively limited 
exposure, which is proportional to the detector surface 
area and experiment lifetime. 


Here we have been able to extend existing constraints 
on nuclearites [31] . based on resonant bar gravitational 
wave detectors, to place new constraints on macro dark 
matter coupling to baryons. If Macros also couple to 
photons with their geometric cross section then the cos¬ 
mic microwave background (CMB) constraints indicating 
o’x/ATx ^ 4.5 X 10“^ cm^/g (dashed line in Fig. [^, when 
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